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'■ Abstract 

Simple multistage closed- (box+reservoir) (MCBR) models of chemical 
evolution, formulated in an earlier attempt, are extended to the limit 
of dominant gas inflow or outflow with respect to gas locked up into 
long-lived stars and remnants. For an assigned empirical differential 
oxygen abundance distribution (EDOD), which can be linearly fit- 
ted, a family of theoretical differential oxygen abundance distribution 
(TDOD) curves is built up with the following prescriptions: (i) the ini- 
tial and the ending point of the linear fit are in common to all curves; 
(ii) the flow parameter, k, ranges from an extremum point to ^foo, 
where negative and positive k correspond to inflow and outflow, re- 
spectively; (iii) the cut parameter, £o, ranges from an extremum point 
(which cannot be negative) to the limit, (Co)oo, related to |k| — > +00. 
For curves with increasing ^o, the gas mass fraction locked up into 
long-lived stars and remnants is found to attain a maximum and then 
decrease towards zero as \k\ — > +00, while the remaining parameters 
show a monotonic trend. The theoretical integral oxygen abundance 
distribution (TIOD) is also expressed. An application is performed 
to the EDOD deduced from two different samples of disk stars, for 
both the thin and the thick disk. The constraints on formation and 
evolution are discussed in the light of the model. The evolution is 
tentatively subdivided into four stages, namely: assembling (A), for- 
mation (F), contraction (C), equilibrium (E). The EDOD related to 
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any stage is fitted by all curves where < Co < (Co)oo for inflowing 
gas and (Co)oo < Co < 1-2 for outflowing gas, with a single exception 
related to the thin disk (A stage), where the range of fitting curves 
is restricted to 0.35 < Co < (Co)oo- The F stage may safely be de- 
scribed by a steady inflow regime (k = —1), implying a flat TDOD, 
in agreement with the results of hydro dynamical simulations. Finally, 
(1) the change of fractional mass due to the extension of the linear 
fit to the EDOD, towards both the (undetected) low-metallicity and 
high-metallicity tail, is evaluated and (2) the idea of a thick disk-thin 
disk collapse is discussed, in the light of the model. 
keywords - galaxies: evolution - stars: formation; evolution. 

1 Introduction 

The empirical metallicity distributions of long-lived stars, belonging to dif- 
ferent populations, constrain models for the formation and the evolution of 
the Galaxy. Simple multistage closed- (box+reservoir) (MCBR) models for- 
mulated in an earlier attempt (Caimmi 2011a, hereafter quoted as Cll) make 
a useful tool in the description of galactic chemical evolution. In particular, 
MCBR models allow for gas inflow or outflow with specified metal abundance 
ratio with respect to the pre existing gas, but general conclusions can also be 
inferred (Cll). The special case of steady inflow regime, where the infalling 
gas exactly balances the net amount of gas turned into stars, finds a counter- 
part in the results of hydrodynamical simulations, where quasi equilibrium is 
attained between inflowing gas, outflowing gas and gas lost via star formation 
(e.g., Finlator and Dave 2008; Dave et al. 2011a,b, 2012). 

For different populations within the Galaxy, the empirical differential 
oxygen abundance distribution (EDOD) shows the existence of three different 
regimes of chemical evolution, where gas inflow is initially the dominant 
process, followed by a nearly steady inflow regime, while gas outflow prevails 
at late times (Cll and further references therein). Understanding different 
stages of Galactic (or sub Galactic) chemical evolution can be used to get 
further insight on the evolution of important global properties of stellar and 
gaseous component of galaxies. 

The current paper is aimed to (i) extension and improvement of MCBR 
models and (ii) application to the thin and the thick disk. The extension of 
the model is concerned with the limit of dominant inflowing or outflowing gas 
rate in comparison to star formation rate and calculation of related quantities 
which, for assigned input values, restricts the parameter space. In addition, 
the theoretical integral oxygen abundance distribution (TIOD) is explicitly 
expressed. The improvement of the model consists in a different choice of 
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the family of fitting theoretical differential oxygen abundance distribution 
(TDOD) curves with respect to an earlier attempt (Cll), which can be used 
for any kind of environment. 

With regard to an assigned stage of evolution, the starting and the ending 
point of the linear fit to the EDOD are in common to each curve of the family, 
instead of the intercept and the ending point as in the parent paper (Cll). 
The difference is negligible if the model is applied to populations with initial 
oxygen abundance close to zero such as the halo (Cll), while it is consistent 
for populations with initial oxygen abundance substantially larger than zero 
such as the thin and the thick disk, which are dealt with in the current paper. 

The EDOD for the thin and the thick disk is inferred from two different 
samples in Section 2. Extension and improvement of MCBR models is per- 
formed in Section 3. An application to the thin and the thick disk is shown 
in Section 4. The results are discussed in Section 5. The conclusion is drawn 
in Section 6. 



2 The EDOD for the thin and the thick disk 



2.1 General remarks 

The abundance distribution of a generic nuclide, Q, for star samples, is binned 
in [Q/H]= log(nq/nH)— log(nq/nH)G, where tiq, nu, are number abundances. 
On the other hand, the comparison between the empirical differential metal 
abundance distribution and its theoretical counterpart predicted by a model 
implies the knowledge of the normalized mass abundance, 0q = Zq/(Zq) . 
It can be seen that log0Q =[Q/H] to a good extent. For further details refer 
to an earlier attempt (Caimmi 2007). 

In terms of mass abundance, the bin centre and the bin semiamplitude 
read: 



which implies a variable bin width in 0q for a constant bin width in [Q/H]. 

The differential metal abundance distribution has been dealt with in ear- 
lier attempts in both non normalized (Pagel 1989; Malinie et al. 1993) and 
normalized (Rocha-Pinto and Maciel 1996; Caimmi 2000, 2001a,b, 2007) 
form, the last one to be used in the following. The related bin centre and 



0Q = T^Pio [Q/H] + ex Pio [Q/h] } ; 
A^ Q = i{exp 10 [Q/H] + - exp 10 [Q/H]-} ; 



(2) 



(1) 
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the bin semiamplitude read: 



V>Q = log 



AN 



(3) 



NA<p Q 



(AN) 



1/21 



A^Vq = log 1 =F 



AN 



(4) 



where AiV is the number of sample stars within the metal abundance bin, 
A(f>Q = A~(()q + A + (f>Q, centered on <pq, N is the total number of sample 
stars, and the uncertainty on AiV has been evaluated from Poissonian errors 
as oajv = (AN) 1 / 2 (e.g., Ryan and Norris 1991). It is worth noticing the 
semiamplitudes, A t ^q, are different and, in particular, A~^q — > — oo for 
AiV = 1, which implies caution has to be used with bins containing a single 
star. For further details refer to earlier attempts (Caimmi 2001b, 2007). 

From this point on, attention shall be restricted only to oxygen, Q=0, 
where 0o — <P an d V'o — ip shall be used to simplify the notation. On the 
other hand, the results of the current paper hold for any primary element 
which is mainly synthesized within type II supernova progenitors. 



Oxygen abundance determination in stellar atmospheres is intrinsically dif- 
ficult, in that different methods yield different results (e.g., Ramirez et al. 
2007; Fabbian et al. 2009) and no general consensus still exists. For this rea- 
son, the population of available samples does not exceed a few hundredths 
at most. For further details refer to earlier attempts (Caimmi 2007, 2008; 
Caimmi and Milanese 2009; Cll). 

With regard to the disk, the data shall be inferred from samples studied 
in earlier attempts (Ramirez et al. 2007; Petigura and Marcy 2011), where 
oxygen abundance has been determined both in presence and in absence 
of the local thermodynamical approximation (LTE) in the former case and 
only in presence of LTE in the latter. The samples extracted from the above 
mentioned parent papers are listed in Table [TJ 

With regard to Ra07 samples, the EDOD has been inferred in an earlier 
attempt (Caimmi and Milanese 2009) both in presence and in absence of 
the LTE approximation and only the latter case shall be considered in the 
current investigation. 

The EDOD inferred from PM11 samples is presented in Table [21 



2.2 The data 
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Table 1: Samples used for determining the empirical differential oxygen abun- 
dance distribution (EDOD) of the thin and the thick disk. Sample denom- 
ination relates to the parent sample (Ra07 - Ramirez et al. 2007; PM11 - 
Petigura and Marcy 2011) and to the population [k - thick disk (TK); n - 
thin disk (TN); stars uncertain if TK or TN are denoted as UN]. Presence 
or absence of the LTE approximation for oxygen abundance determination 
is denoted as yes or no, respectively. The total number of sample stars is N. 



sample 


population 


N 


LTE 


Ra07k 


TK 


133 


no 


Ra07n 


TN 


310 


no 


PMllk 


TK + UN 


14 + 20 


yes 


PMlln 


TN 


624 


yes 



2.3 Linear fit to the EDOD 

Both Ra07 and PM11 samples are biased towards low metallicities, [Fe/H] ~ 
— 1 which, using the [0/H]-[Fe/H] relation inferred in the parent paper 
(Ramirez et al. 2007), translates into ~ 0.3 for both the thin and the 
thick disk. On the other hand, thin disk stars exhibit oxygen abundance 
above 0.3 within both Ra07n (Caimmi and Milanese 2009) and PMlln (Ta- 
ble^ samples, while both Ra07k (Caimmi and Milanese 2009) and PMllk 
(Table [2]) samples show a low-metallicity tail below 0.3. 

It may safely be assumed the low-metallicity tail relates to a negligible 
mass fraction for the thin and (to a lesser extent) the thick disk. Accordingly, 
sufficiently large samples listed in Table [T] (PMlln, Ra07n, Ra07k) can be 
considered as representative for the thin or the thick disk, while the remaining 
low sample (PMllk) is added for comparison only. 

For all samples listed in Tabled! the EDOD is plotted in Fig.[lJ where up- 
per and lower panels represent PM11 and Ra07 samples, respectively. Points 
with lower error bars attaining the horizontal axis (including the one out of 
scale) correspond to bins containing a single star. 

The main feature of the EDOD, plotted in Fig.[H is the presence of up 
to four regions exhibiting a nearly linear trend, which shall be named as 
A, F, C, E, for increasing oxygen abundance. Linear fits to each region, 
ip = acf) + b, are performed using the bisector regression (Caimmi 2011b), 
leaving aside bins related to single stars, unless the total number of points 
reduces to two. Aiming to a unified description, the number of regions shall 
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Table 2: The empirical, differential oxygen abundance distribution (EDOD) 
in the thin and thick disk, deduced from the PM11 sample (N = 658). The 
error on the generic bin height has been estimated from the Poissonian error. 
Disk stars belonging to uncertain population (N = 20) have been included 
within thick disk stars (N = 14) to get larger population (N = 34). The bin 
width in [0/H]= log0 is ^0.05. Label captions: TN - thin disk; TK - thick 
disk; UN - uncertain if belonging to thin or thick disk. See text for further 
details. 





TN 


TK + UN 


AN 














TN 


TK 


UN 


1.2673D- 


■1 






+5.3464D- 


3 





1 





1.5954D- 


1 



















2.0085D- 


1 



















2.5286D- 


1 



















3.1833D- 


1 



















4.0075D- 


1 


-1.4573D- 





-4.9465D- 


1 


2 


1 





5.0451D- 


•1 


-1.0802D- 


■0 


-2.9362D- 


1 


6 


1 


1 


6.3514D- 


•1 


-4.9596D- 


1 


-3.9362D- 


1 


29 


1 


1 


7.9960D- 


1 


-3.7170D- 


■2 


-1.9259D- 


1 


105 


2 


2 


1.0066D- 


■0 


-1.8480D- 


-2 


-2.9259D- 


1 


138 


1 


3 


1.2673D- 





-7.9393D- 


-3 


-9.1564D- 


2 


178 


2 


6 


1.5954D- 





-2.7200D- 


-1 


-1.9156D- 


1 


122 


5 


3 


2.0085D- 


-0 


-8.5630D- 


1 


-7.1753D- 


1 


40 





3 


2.5286D- 





-1.9563D- 


-0 


-1.2947D- 





4 





1 
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Table 3: Regression line slope and intercept estimators, a and b, and related 
dispersion estimators, <7a, and ay for bisector regression models applied to the 
oxygen abundance distribution (EDOD) plotted in different panels of Fig.[TJ 
The method has dealt with each region (X) separately. Data points on the 
boundary between adjacent regions are used for determining regression lines 
within both of them. 



Y 
A 




a 



























sample 


A 


+3. 


.6420 


E+0 


1 


.4136 


E- 


-1 


-2 


.8981 


E+0 


6 


.6652 


E- 


-2 


PMlln 


F 


+6. 


.3321 


E-2 


5. 


.7541 


E- 


-3 


-8 


.6073 


E-2 


5 


.8764 


E- 


-3 




C 


-1 


.1689 


E+0 


7. 


.0570 


E- 


-2 


+1. 


.5193 


E+0 


1 


.2997 


E- 


-1 




E 


-1 


.8294 


E+0 


8. 


0990 


E- 


-2 


+2. 


.7116 


E+0 


1 


.8562 


E- 


-1 




A 


+2 


.1371 


E+0 


3. 


9906 


E- 


-2 


-1 


.2771 


E+0 


2 


.2528 


E- 


-2 


Ra07n 


F 


-4. 


.3676 


E-l 


3. 


.0185 


E- 


-2 


+3. 


.4846 


E-l 


2 


.8952 


E- 


-2 




CE 


-1 


.0703 


E+0 


1. 


.3474 


E- 


-2 


+1. 


1649 


E+0 


2 


.2887 


E- 


-2 




AF 


+2. 


.7877 


E-l 


7. 


.8677 


E- 


-2 


-5 


.1247 


E-l 


7 


.1999 


E- 


-2 


PMllk 


CE 


-1 


.1799 


E+0 


1. 


.8906 


E- 


-2 


+1. 


.6773 


E+0 


4 


.1322 


E- 


-2 




A 


+1. 


.8383 


E+0 


2. 


.4630 


E- 


-1 


-1 


.1545 


E+0 


1 


.3187 


E- 


-1 


Ra07k 


F 


-6 


.4249 


E-2 


2. 


.8011 


E- 


-3 


+8. 


.2504 


E-3 


2 


.3791 


E- 


-3 




CE 


-8 


.1530 


E-l 


3. 


.4051 


E- 


-2 


+7. 


.8279 


E-l 


4 


.4835 


E- 


-2 





be kept equal to four in all cases and adjacent regions with coinciding linear 
trend shall be indicated by related letters e.g., CE means regions C and E 
are fitted by the same line. The regression line slope and intercept estimators 
and related dispersion estimators are listed in Table [3] for each region of the 
EDOD plotted in Fig.HJ 

The regression lines are represented in Fig. [2] for each region, according 
to slope and intercept values listed in Table [3J To ensure continuity, ad- 
jacent regions are bounded by the intersections of related regression lines. 
The results are listed in Table HI where O denotes oxygen abundance ranges 
without data, < < <pi and > 0/, where <pi and <pf are the minimum and 
maximum oxygen abundance, respectively, within sample stars. Accordingly, 
a vertical line instead of a regression line is considered for the intersection 
points related to O-A and E-0 transitions. 

In conclusion, the EDOD related to the thin and the thick disk may be 
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Table 4: Transition points between adjacent regions, as determined from the 
intersection of related regression lines, for the oxygen abundance distribution 
(EDOD) plotted in different panels of Fig.[TJ The repetition of the data 
passing from a stage to the next one means that the interpolation line remains 
unchanged and any point may be considered as a transition point. Cases on 
the left and on the right relate to thin and thick disk, respectively (left and 
right panels of Fig. [I]). Cases at the top and at the bottom relate to PM11 
and Ra07 samples, respectively (top and bottom panels of Fig. [I]). 



sample 


PMlln / Ra07n 


PMllk / Ra07k 


transition 


i> 


i> 


O-A 


4.0075 E-l 


-1.4386 


E+0 


5.0451 E- 


-1 


-3.7183 


E- 


1 


A-F 


7.8579 E-l 


-3.6316 


E-2 


5.0451 E- 


-1 


-3.7183 


E- 


1 


F-C 


1.3028 E+0 


-3.5782 


E-3 


1.5012 E- 


-0 


-9.3977 


E- 


2 


C-E 


1.8052 E+0 


-5.9085 


E-l 


1.5012 E- 


-0 


-9.3977 


E- 


2 


E-0 


2.5286 E+0 


-1.9141 


E-0 


2.5286 E- 


-0 


-1.3061 


E- 





O-A 


4.0075 E-l 


-4.2062 


E-l 


2.5286 E- 


-1 


-6.8969 


E- 


1 


A-F 


6.3155 E-l 


+7.2630 


E-2 


6.1116 E- 


-1 


-3.1016 


E- 


2 


F-C 


1.2887 E-0 


-2.1439 


E-l 


1.0313 E- 


-0 


-5.8008 


E- 


2 


C-E 


1.2887 E-0 


-2.1439 


E-l 


1.0313 E- 


-0 


-5.8008 


E- 


2 


E-0 


2.5286 E-0 


-1.5414 


E-0 


2.0085 E- 


-0 


-8.5474 


E- 


1 



8 



conceived, to a satisfactory extent, as due to the contribution of four (A, F, 
C, E) or three (A, F, CE) or two (AF, CE) regions within which the trend is 
linear. An interpolation in terms of MCBR models of chemical evolution is 
highly attractive in that the corresponding TDOD shows, as a special case, 
a linear trend (Cll). 



3.1 General remarks 

Simple MCBR models have been presented in an earlier attempt (Cll) and 
an interested reader is addressed therein for an exhaustive formulation. Only 
what is relevant for the application shall be repeated here and related ex- 
tension and improvement shall be performed. The main assumptions of the 
model are listed below. 

• System structured as a box and a reservoir where gas can be exchanged 
between the two but mass conservation holds. 

• Instantaneous recycling within the box, where stars are divided into 
two categories, namely (a) short-lived, which instantaneously evolve, 
and (b) long-lived, all of which are still evolving. 

• Instantaneous mixing within the box. 

• Gas outflow from the box into the reservoir or inflow into the box from 
the reservoir, at a rate proportional to the star formation rate and with 
fixed composition. 

• Inhibition of star formation within the reservoir. 

• Oxygen abundance proportional to metal abundance, according to the 
forthcoming Eq. (jSj). 

In this picture, the gas mass fraction, fi, (normalized to the initial mass) 
and the oxygen abundance, 0, (normalized to the solar value) are related as: 



3 The model 




n (Z o ) /(cp") 



(5) 



C = 



P 



[A p - k(1 - Co)] 



(6) 



Z = A Z 



(7) 
(8) 
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where (Zq) q is the solar oxygen abundance, p is the yield, p" the effective 
yield, k the flow (positive for outflow and negative for inflow), Co the cut, 
Zo the oxygen abundance, Z the metal abundance, Aq> = 2 to an acceptable 
extent and the index, i, denotes values at the beginning of evolution. 

More specifically, k represents the ratio of gas mass locked into long-lived 
stars and stellar remnants to gas mass outflowed into the reservoir or inflowed 
from the reservoir and Co represents the ratio of oxygen abundance within 
the flowing gas to oxygen abundance within the pre existing gas. For further 
details refer to the parent paper (Cll). In principle, — oo < k < +00 and 
< Co < Ao/(Zo)f. The domain, Co > Ao/(Zo)f, is in contradiction with 
the assumption of fixed Co during the evolution. 

The TDOD predicted by the model reads: 



(Zoh 



- 1 



cp" 



b = log 



loefl 



'Oj© 



+ b 



(Zoh 



cp" 



log(l - C0i 



(9) 
(10) 



where b is the intercept of the curve on the {0<j)ip) plane. For further details 
refer to the parent paper (Cll). 



3.2 Theoretical integral oxygen abundance distribu- 
tion (TIOD) 

An explicit expression of the TIOD could be needed for a number of appli- 
cations. The theoretical counterpart of the EDOD, via Eq. (E]), reads: 



ij) = log 



dN 



(N f - Ni) d0 



(11) 



where Ni, Nf, denote the total number of long-lived stars with normalized 



oxygen abundance, < 



< 



respectively. 



The combination of Eqs. (jH]) and fllip yields after some algebra: 
diV 



10TL-cc 



v(Z o ) /(cp")-l 



N f -N 

which can be integrated. The result is: 



d0 



N-Nj 
N f -N 



-10 



b P" 



(Z 



o)g 



v(Z o ) /(cp") 



:i-cc 



v(Z o ) /( C p") 



(12) 



(13) 



where N is the total number of stars with normalized oxygen abundance not 
exceeding 0. Accordingly, the TIOD is expressed by Eq. ffTBl . 
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3.3 A family of TDOD curves 



The TDOD, expressed by Eq. ©, via ©, ©, ©, (HDD, depends on the in- 
put parameters, [Ao,p, (Zo)o, 0«, 0/, /ii] and the output parameters, (k, Co)- 
With regard to a selected stage of evolution (or adjacent stages fitted by the 
same line), let (</>j,^>j) and ((f>f,i/jf) be the starting and the ending point, 
respectively, of the related linear fit. The family of TDOD curves, passing 
through the above mentioned points, necessarily satisfy the relations: 



1pf 



oj© 



- 1 



cp 

[Zoh 
cp" 



- 1 



log(l - C0 4 ) + b 
log(l — o 



+ b 



(14) 
(15) 



which makes a system of two equations in the two unknowns, (k, (q) 
The combination of Eqs. (1141) and (1151) yields after some algebra: 



cp 



(Z 



o)o 



l + (V»/--0i) |l0g- 



\ — cc 



(16) 



which, for assigned Co, can be solved by repeated iterations using Eq. ([6]) 
until — k^ 1 ^ < e, where n is the last iteration and e a fixed threshold. 
Then the generic curve of the family is defined by the parameters, (k, Co)- 
The intercept, b, is determined via Eq. (fT4"j) or ffTSI) . 

Both the EDOD and the TDOD can be represented on the {0<pip) plane, 
the former and the latter normalized to the long-lived star population of the 
sample and the environment under consideration, respectively, which trans- 
lates into a normalization constant, logCN- More specifically, the TDOD has 
to be vertically shifted on the (Ocpip) plane by a value, logCN, for matching 
to the EDOD. The result is: 



10' 



v m - 



o)q 



,(Z O ) /( C p") 



(17) 



where b is defined by Eq. (flQj) . For further details refer to the parent paper 
(Cll), keeping in mind the family of TDOD curves defined therein have 
in common the points, (0,6) and (0/,^/), instead of (4>i,ipi) and (4>f,ipf), 
respectively. 



3.4 The linear limit 

In the special case, c = 0, Eq. 



c« 



o 



reduces to: 
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which implies Co > 1 f° r K < 0; < Co < 1 f° r K > ^4oP; Co — 1 f° r 
|«| — > +oo; Co < for < k < A p. The last subdomain has to be excluded 
as negative cut values, by definition, have no physical meaning. 

With regard to the remaining domain, Eqs. ©, ([9]), (|T0i) . (|T3i) . reduce to: 



— = exp 

Hi 

ip = a(fi + b ; 

fit 



p" 



log 



1 {Zo)e 
a = 

In 10 p" 
N — Ni 11 



lnl0)a 



a(j>i 



[exp lo (a0 + b) - exp lo (a0j + b)} 



(19) 
(20) 
(21) 

(22) 
(23) 



N f - Ni In 10 a 

implying a linear TDOD on the {0<j)ip) plane. Accordingly, the limit, c — > 0, 
shall be quoted as the linear limit. 

Keeping in mind the starting and the ending point of the linear fit are 
(0i, ■?/>«), {4>fi i>f), respectively, an empirical counterpart of Eq. (120]) reads: 



if) = ipi + 

<Pf ~ <Pi 

which implies, by comparison with Eq. f l2"U|) via ([7]) and 



cj) f 



a 



a (l + k) 



In 10 p 

where ao is the slope related to stagnation regime, k 
n = a /do — 1 and Eq. (fl8|) takes the expression: 

Co - 1 + -, / — , 

1 — a/ ao 

where the output parameters, (k, Co), ma Y be determined using Eqs. (121)]) 



(24) 

(25) 
(26) 

0. Accordingly, 
(27) 



3.5 The dominant flow limit 



In the limit of negligible star formation rate with respect to flow rate, \k\ — > 
+oo, the following relations hold after some algebra: 



lim 



cp" 



(2f 



12 



lim 

\k\— >+oo 

lim 
lim 

\k\— >+co 



j O)& 



p" 1 — O 

1 — cc 



I -Co 



log 
64 



1 - C(pi_ 

(Zoh 



log 



cp 



log(l - C6 



log 



Mi 



1 1 



ja - n f 1 - Co ■ 



(29) 
(30) 
(31) 



and Eqs. © and Q) reduce to: 



\k\-*+oo fa \(pij 

lim ip = — log — + log 



i-C 



o 



Hi 



Hi - fJL f 1 - Co 0i. 



(32) 
(33) 



where the last relation implies a divergence for the intercept, i/)(0) — > oo. 
The particularization of Eq. (I3"3"j) to the points, (<^,^), (0/, ■?/>/), yields: 



/ / Co -2 , 
iPf-n = ~ x — log 



c 



o 



accordingly, the cut, Co ; reads after some algebra: 

(V>/ - V.) + 2 log(0//0i) _ a/ 7 + 2 



C 



o 



7 



- Vi) + log(0//0i) a/7 4- 1 
log(0//0i) 



(34) 

(35) 
(36) 



where a is the slope of the linear fit, expressed by Eq. f[2"5"j) . and 7 > 0. 
The intercept, 6, according to Eq. (13TI) . can be expressed as: 



6 = 6'- 
6' = log 



(Zoh 
cp" 

Hi 



l0g(l - C0j) 

I 1 



Mi - m/ 1 - c 



o <Pi 



(37) 
(38) 



where k — >■ =Foo necessarily implies Hi > an d) i n turn, Co < 1 5 within the 
domain of the reduced intercept, b', expressed by Eq. (138]) . Keeping in mind 
Eq. ([6]), the last term on the right-hand side of Eq. (137|) in the case under 
discussion shows a divergence as: 



lim 

\k\— >+oo 



(Zoh 1 



Cp 



log(l - c0i) 1 = fr* 2 lim log(l - c0j) = 00 ; 

I 1 — (,0 |«|-H-oo 



(39) 
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accordingly, the normalization constant, C*n, must be related to the reduced 
intercept, b', whose empirical counterpart equals ipi via Eqs. (fUj) and (1371) . 
Then Eq. (|T7|) translates into: 



Hi 

which is equivalent to ipi = b' + logCn. 



c N = 10*^—^(1 -Co)0i ; (40) 



4 The (O< ) plane 



The linear limit, c = 0, is represented on the (O/tCo) plane, via Eq. ffTSl) . as 
an equilateral hyperbola of asymptotes, k = and Co — 1> shown in Fig. [31 
It can be seen c > for points lying between the branches of the hyperbola, 
while c < for points lying outside either branch of the hyperbola. 

The steady inflow regime, k = —1, is represented on the (OkCo) plane as 
a vertical axis. It can be seen fif > fa in the strong inflow regime, k < — 1, 
and (if < fa in the weak inflow regime, — 1 < k < 0, or in the stagnation 
regime, k = 0, or in outflow regime, k > 0, while /!/ = /ij in the steady inflow 
regime. 

The limit of complete gas exhaustion, via Eq. ([5]), reads c = 1/0/ which, 
using Eq. (jH]), can be cast under the explicit form: 



o 



A p- 



j o)g 



1 

K 



(41) 



that is an equilateral hyperbola with asymptotes, k — 0, £o — 1> as shown 
in Fig.[3j It can be seen /!/ > for points lying between the branches of 
the hyperbola, while /i/ < for points lying outside either branch of the 
hyperbola. The latter implies lack of physical meaning, together with points 
lying on the negative semiplane, Co < 0. Accordingly, the parameter space 
on the (O/tCo) plane is restricted to the region bounded by the branches of 
the hyperbola, expressed by Eq. fj4T]) . and the horizontal axis, Co = 0, where 
> 0, Co > 0, as shown in Fig.[3j 

It is worth noticing c0/ — > 1 implies \i/j((j)f) \ — > +oo via Eq. ([9]) and the 
related curve cannot belong to the family, defined by Eqs. f[T4"|) and fTTS]) . 

5 Results 

The EDOD related to the samples considered in Section 2 can be divided 
into four, three, or two regions, where the trend is linear to a good extent, 
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as shown in Fig.[2J In the light of the model, each region corresponds to a 
different stage of chemical evolution, which can tentatively be related to a 
different stage of dynamical evolution. For further details refer to the parent 
paper (Cll). 

The fractional mass of the box (with respect to the initial value) attains 
the maximum value at the end of the last stage, P, where the TDOD still 
shows a slope, a > a , i.e. inflow regime (into the box from the reservoir). 
The related value is (/xp)/ + (sp)/- The fractional stellar mass of the box 
at the end of evolution, which coincides with the end of the last stage, L, is 
(sl)/. Accordingly, the mass ratio of the box at the end of evolution to the 
outflowed gas reads: 

M box = (/j L )/ + (s L )f . (42) 

M ofl (ji P ) f + (sp) f - (ji L ) f - (s L ) f ' 

provided earlier stages with respect to P are in inflow regime and later stages 
in outflow regime. 

Following the procedure outlined in the parent paper (Cll) and keeping 
the same values of the input parameters, (Z ) Q = 0.0056, p/ (Z ) Q = 1.0340, 
Hi = 1, Si = 0; Di = 0, the normalization constant, the cut, the flow, the 
active gas mass fraction, the star mass fraction and the inflowed or outflowed 
gas mass fraction, at the end of each stage of evolution, can be computed. 
The results are listed in Tables |5] and El where different cases are related 
to different EDOD data plotted in Fig.[TJ The index, U, denotes a stage 
(or adjacent stages) of evolution where the trend is linear to a good extent 
and the indexes, lin and oo, mark the linear and the dominant flow limit, 
respectively. 

The mass ratio of the box at the end of evolution to the outflowed gas, 
expressed by Eq. (1421) . and the box mass fraction at the end of evolution, 
(M box )j/(M box )j = (//l)/ + ( s l)/, can be inferred from Table El where P = 
F, A, AF, and L = E, CE, according to the case. The result is (leaving aside 
the poorly populated PMllk sample) Mb ox /M fl =1.4-2.6 for the thin disk 
and 2.9 for the thick disk; (Mb ox )//(^box)i =2.5-25.0 for the thin disk and 
5.6 for the thick disk. 

The validity of the linear limit implies cut values very close to unity, 
Co ~ 1, unless k — > 0" or k — » (A p) + according to Eq. (fTSl) . as can be seen 
in Table El On the other hand, inflowing gas from intergalactic medium is 
expected to be oxygen-poor with respect to galaxies, while outflowing gas 
into galactic medium is expected to be slightly oxygen-rich with respect to 
pre existing gas. For this reason, the departure from the linear limit must 
also be considered. 
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Table 5: Output parameters, [(Co)u]im and [(Co)u]oo, related to the linear 
limit (lin) and to the dominant flow limit (oo), for simple MCBR models 
where the theoretical differential oxygen abundance distribution (TDOD) 
provides a linear fit to the empirical differential oxygen abundance distribu- 
tion (EDOD) plotted in different panels of Fig.[TJ Four, three, or two stages 
of evolution are considered, according to the linear trends exhibited by the 
EDOD. 



u 


[(Co)u]lin 


[(Co)u]oo 


sample 


A 


1.0012E+0 


1.1726E+0 


PMlln 


F 


1.0101E+0 


1.8702E+0 




C 


9.9351E-1 


6.8213E-1 




E 


9.9655E-1 


8.7566E-1 




A 


1.0019E+0 


1.2860E+0 


Ra07n 


F 


7.0971E-1 


1.4628E+1 




CE 


9.9252E-1 


7.1698E-1 




AF 


1.0070E+0 


1.6302E+0 


PMllk 


CE 


9.9360E-1 


7.7029E-1 




A 


1.0023E+0 


1.3778E+0 


Ra07k 


F 


1.0138E+0 


2.1469E+0 




CE 


9.8770E-1 


4.2840E-1 
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Table 6: Output parameters, (C\j)n, k u> (/iu)/> ( s u)/, (-Du)/) f° r simple 
MCBR models where the theoretical differential oxygen abundance distri- 
bution (TDOD) provides a linear fit to the empirical differential oxygen 
abundance distribution (EDOD) plotted in different panels of Fig. [I] Other 
captions as in Table [5j 

U (Cu)n «u (a*u)/ (sy)f (Dv)f 



A 
F 
C 
E 



1.0534E- 
4.9391E- 
2.7317E- 
5.8008E- 



-9.6714E+0 
-1.1508E+0 
+1.7831E+0 
+3.3557E+0 



2.5253E+1 
2.7230E+1 
7.0432E+0 
3.3456E-1 



2.7969E+0 
1.5911E+1 
2.3164E+1 
2.4704E+1 



-2.7050E+1 
-4.2141E+1 
-2.9207E+1 
-2.4039E+1 



A 1.6306E-1 
F 5.6841E-1 
CE 2.3602E-1 



-6.0883E+0 
+3.9895E-2 
+ 1.5483E+0 



3.1135E+0 
1.6078E+0 
7.5727E-2 



4.1536E-1 
1.8633E+0 
2.4645E+0 



-2.5289E+0 
-2.4711E+0 
-1.5402E+0 



AF 5.9299E-1 -1.6637E+0 1.8961E+0 1.3500E+0 -2.2461E+0 
CE 8.7126E-1 +1.8093E+0 1.1633E-1 1.9835E+0 -1.0999E+0 



A 1.7169E-1 
F 5.5095E-1 
CE 9.4260E-1 



-5.3769E+0 
-8.4703E-1 
+9.4119E-1 



4.5570E+0 
4.2824E+0 
6.8384E-1 



8.1266E-1 
2.6077E+0 
4.4615E+0 



-4.3696E+0 
-5.8901E+0 
-4.1454E+0 
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The results are shown in Fig.H]where, in addition to the linear limit shown 
in Fig. 121 the following cases are also represented: Co = 0, (Co)oo, for k < i.e. 
inflow regime and Co — (Co)ooj 1-2, for k > i.e. outflow regime. In general, 
curves from up to down relate to decreasing cut, Co- All curves maintain 
close to the related linear limit and fit to the data, with the exception of the 
A stage corresponding to the PMlln sample, where curves within the range, 
< Co < 0.35, do not match the whole set of error boxes. The curve related 
to Co = 0.35 is also plotted on the top left panel of Fig.Hl 

With regard to the PMlln sample restricted to the A stage, values of 
parameters which characterize the family of TDOD curves, considered in 
Subsection 3.2, are listed in Tabled The linear limit, the dominant flow limit, 
together with other curves plotted in Fig.[U top left panel, A stage, are also 
included. Due to the occurrence of inflow regime, the flow, k, and related 
fractional flowing mass, D, are negative by definition. All the parameters 
listed in Table [7] exhibit a monotonic trend with the exception of the long- 
lived star fractional mass, Sf, which attains a maximum at Co ~ 0.70 and 
decreases to zero in the dominant flow limit, Co — (Co)oo- in particular, the 
TDOD fits to the EDOD within the range, 0.35 ~ Co < (Co)oo, related to 
— 1.85 ~ k > — oo, 3.77 ~ fif ~ 49.52, < Sf ~ 3.63. It can also be seen 
the normalization constant, Cn, is changed only slightly with respect to the 
linear limit, Co = (Co)iin- 

In the general case of a varying oxygen abundance within the flowing gas 
i.e. varying cut, (Co)min < Co < (Co) max , the TDOD belonging to the family 
of curves passing through the points, (4>i,ipi), (<pf,ipf), intersects all curves 
characterized by constant Co i n the range under consideration, lying within 
the region bounded by the curves where Co = (Co)min and Co = (Co)max- As 
shown in Fig.Hl the TDOD fits to the EDOD over the whole range of Co> 
with the restriction, (Co)min = 0.35, for the A stage, top left panel (middle 
full curve). 

Values of active gas, long-lived star and flowing gas mass fraction at 
the end of each stage are listed in Table [8] for the linear limit (upper lines) 
and the lowest curve fitting to the data (lower lines) which corresponds to 
(Co)min = 0.35 for the PMlln sample and (Co)min = for Ra07n, PMllk, 
Ra07k samples, with regard to the initial A stage. The linear limit is assumed 
during the following F, C and E stage. An inspection of Table [H] discloses 
that the active gas mass fraction is lowered by a factor up to about 7, while 
long-lived star and flowing gas mass fraction are lowered by a factor up to 
about 4, passing from the linear limit, Co = (Co)im, to the lowest curve fitting 
to the data, Co — (Co)min, at the end of evolution. 
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Table 7: Values of the cut, Co, the flow, k, the fractional active gas mass, 
fif, the fractional long-lived star mass, Sf, the fractional flowing gas mass, 
Df, the parameter, c, and the normalization constant, Cn, for the family 
of theoretical differential oxygen abundance distributions which contains the 
curves plotted in Fig.lH top left panel restricted to A stage. In addition, 
Ciin = 1.0012 and Coo = 1.1726, as listed in Table [5j For further details refer 
to the text. 



Co 


— K 




s f 


-Df 


c 


Cn 




0.00 


1.2587E+0 


1.7268E+0 


2.8099E+0 


3.5367E+0 


+1.2284E+0 


5.3730E- 


-2 


0.10 


1.3860E+0 


2.1334E+0 


2.9360E+0 


4.0694E+0 


+ 1.2176E+0 


5.6624E- 


-2 


0.20 


1.5407E+0 


2.6596E+0 


3.0693E+0 


4.7290E+0 


+1.2032E+0 


5.9859E- 


-2 


0.30 


1.7323E+0 


3.3487E+0 


3.2072E+0 


5.5560E+0 


+1.1839E+0 


6.3482E- 


-2 


0.35 


1.8462E+0 


3.7728E+0 


3.2765E+0 


6.0493E+0 


+1.1718E+0 


6.5455E- 


-2 


0.40 


1.9754E+0 


4.2627E+0 


3.3450E+0 


6.6077E+0 


+1.1574E+0 


6.7546E- 


-2 


0.50 


2.2932E+0 


5.4925E+0 


3.4740E+0 


7.9665E+0 


+1.1201E+0 


7.2111E- 


-2 


0.60 


2.7752E+0 


7.1749E+0 


3.5792E+0 


9.7541E+0 


+1.0654E+0 


7.7246E- 


-2 


0.65 


3.0047E+0 


8.2475E+0 


3.6152E+0 


1.0863E+1 


+ 1.0282E+0 


8.0052E- 


-2 


0.70 


3.3446E+0 


9.5204E+0 


3.6341E+0 


1.2154E+1 


+9.8155E-1 


8.3030E- 


-2 


0.75 


3.7663E+0 


1.1040E+1 


3.6292E+0 


1.3669E+1 


+9.2179E-1 


8.6193E- 


-2 


0.80 


4.3029E+0 


1.2839E+1 


3.5919E+0 


1.5456E+1 


+8.4346E-1 


8.9553E- 


-2 


0.90 


5.9747E+0 


1.7756E+1 


3.3683E+0 


2.0125E+1 


+5.8900E-1 


9.6916E- 


-2 


1.00 


9.6022E+0 


2.5143E+1 


2.8067E+0 


2.6950E+1 


+1.1200E-2 


1.0523E- 


-1 


Clin 


9.6714E+0 


2.5253E+1 


2.7969E+0 


2.7050E+1 





1.0534E- 


-1 


1.10 


2.3278E+1 


3.6730E+1 


1.6039E+0 


3.7334E+1 


-2.2400E+0 


1.1463E- 


-1 


1.15 


7.5661E+1 


4.5018E+1 


5.8956E-1 


4.4607E+1 


-1.0964E+1 


1.1977E- 


-1 


1.16 


1.3622E+2 


4.6946E+1 


3.3979E-1 


4.6286E+1 


-2.1067E+1 


1.2084E- 


-1 


1.17 


6.7115E+2 


4.8979E+1 


7.1595E-2 


4.8051E+1 


-1.1033E+2 


1.2192E- 


-1 


Coo 


+oo 


4.9516E+1 





4.8516E+1 


— oo 


1.2219E- 


-1 
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Table 8: Values of active gas, long-lived star and flowing gas mass fractions 
at the end of each stage for the extreme situations of inflowing gas during the 
earlier stage (A) with oxygen abundance equal to Co — (Co)im (upper lines) 
and Co = (Co)min (lower lines) with respect to pre existing gas, in connection 
with the thin (top cases) and thick (bottom cases) disk. The values of the 
cut related to lower curves still consistent with the data, plotted in Fig.Hl are 
(Co)min = 0.35 (PMlln sample), (Ra07n, PMllk, Ra07k samples). During 
the subsequent F, C, and E stage, Co = (Co)im has been assumed. 



u 


Mf 




( s u)/ 


(Aj)/ 


sample 


A 


2.5253E4 


-1 


2.7969E+0 


-2.7050E+1 


PMlln 




3.7728E+ 


-0 


6.2, /DOri/+U 


ex n a nor? i n 




F 


2.7230E+ 


-1 


1.5911E+1 


-4.2141E+1 






4.0681E4 


-0 


5.2358E+0 


-8.3039E+0 




C 


7.0432E4 


-0 


2.3164E+1 


-2.9207E+1 






1.0523E+0 


6.3194E+0 


-6.3717E+0 




E 


3.3456E- 


-1 


2.4704E+1 


-2.4039E+1 






4.9984E- 


-2 


6.5495E+0 


-5.5995E+0 




A 


3.1135E- 


-0 


4.1536E-1 


-2.5289E-0 


Ra07n 




1.2393E- 


-0 


7.7920E-1 


-1.0185E-0 




F 


1.6078E- 


-0 


1.8633E-0 


-2.4711E-0 






6.3995E- 


-1 


1.3555E-0 


-9.9548E-1 




CE 


7.5727E- 


-2 


2.4645E-0 


-1.5402E-0 






3.0142E- 


-2 


1.5948E-0 


-6.2497E-1 




AC 


1.8961E- 


-0 


1.3500E-0 


-2.2461E-0 


PMllk 




3.0081E- 


-1 


1.7505E-0 


-1.0513E-0 




CE 


1.1633E- 


-1 


1.9835E-0 


-1.0999E-0 






1.8456E- 


-2 


1.8510E-0 


-8.6944E-1 




A 


4.5570E- 


-0 


8.1266E-1 


-4.3696E-0 


Ra07k 




1.3223E- 


-0 


1.0925E-0 


-1.4148E-0 




F 


4.2824E- 


-0 


2.6077E-0 


-5.8901E-0 






1.2474E- 


-0 


1.5537E-0 


-1.8011E-0 




CE 


6.8384E- 


-1 


4.4615E-0 


-4.1454E-0 






1.9836E- 


-1 


2.0941E-0 


-1.2925E-0 
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6 Discussion 



As shown in Fig.|2], the EDOD inferred from PM11 and Ra07 samples for 
the thin and the thick disk (see Table [1]) can be fitted by a straight line on 
four, three or two adjacent regions. In the light of simple MCBR models of 
chemical evolution, the TDOD is represented by a family of curves (including 
the straight line) passing through the starting and the ending point, (fa, ipi), 
(4>f,ijjf), respectively, of a selected stage of evolution. Accordingly, different 
stages of evolution relate to different linear fits to the EDOD, where an initial 
strong inflow regime (a > 0) is followed by a steady (a = 0) or moderate 
(—0.42 < a < 0) inflow regime and finally by a stagnation (a = —0.42) or 
outflow (a < —0.42) regime. For further details refer to an earlier attempt 
(Cll). 

The above description is in agreement with the results of hydrodynam- 
ical simulations, where quasi equilibrium is attained between inflowing gas, 
outflowing gas and gas lost via star formation, after an early stage of strong 
inflow and, presumably, before a late stage of strong outflow (e.g., Finlator 
and Dave 2008; Dave et al. 2011a,b; 2012). Within the formalism of simple 
MCBR models, both the inflow and the outflow rate are proportional to the 
star formation rate and, for this reason, only the net effect, k = k-^ + k q_, 
Co = (Co)ifl + (Co)ofi, is considered. 

Though the assumption of null stellar mass fraction holds to an acceptable 
extent for the starting configuration, still the presence of stars with lower 
oxygen abundance than in sample objects, <fi < 0.4 or [O/H] < —0.4 for 
the thin disk and (p < 0.2 or [O/H] < —0.7 for the thick disk, has to be 
considered. To this aim, samples unbiased towards low oxygen abundance 
(not available at present) should be used. 

The thick disk is expected to be affected from this inconvenient to a 
larger extent with respect to the thin disk, due to a more extended low- 
metallicity tail. In any case, the results of the current paper may be thought 
of as representative of the chemical evolution of the thin disk for > 0.4 or 
[O/H] > -0.4 and of the thick disk for > 0.3 or [O/H] > -0.5, the last 
related to [Fe/H] > —1, where the data are unbiased (Ramirez et al. 2007; 
Petigura and Marcy 2011). 

As shown in Fig.Hl the TDOD fits to the EDOD regardless of the value 
of the cut, Co i-e. the oxygen abundance ratio of flowing to pre existing 
gas, within a plausible range where the physical meaning of the problem 
is preserved. The sole exception relates to the PMlln sample (Fig.Hl top 
left), A stage, where curves within the range, < Co < 0.35, provide a bad 
fit to the data, which implies inflow of low-metallicity gas within the thin 
disk could not take place when the pre existing gas was sufficiently oxygen- 
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enriched, [0/H] ~ -0.4. 

The above result can be extended to the case of variable cut, (Co) mm < 
Co < (Co)max, where the TDOD lies between the region bounded by the 
curves related to (Co)min and (Co)max, respectively. Then the differential 
distribution of oxygen abundance in thin disk and thick disk stars (leaving 
aside the low-metallicity tail) depends only slightly on the oxygen abundance 
of the flowing gas. 

Biases towards low-metallicity and high-metallicity stars may be corrected 
under the assumption that the TDOD related to (p < (pi and <fi > <pf, respec- 
tively, allow analytical continuation. Let </> min and max be the unknown 
initial and final normalized oxygen abundance related to stages A and E, 
respectively. 

The fractional number of long-lived stars with normalized oxygen abun- 
dance, m i n < < (pi, can be inferred from Eq. ({TBI as: 



-10° 



Nf-Nt ' (Z Q ) Q 

;i-c A l ) (Zo)e/(cA?i " ) -(l-c A min )( z °^/( CA?5 ")] ; (43) 

or, in the linear limit (c — > 0), from Eq. (123 p as: 
(SN) A 1 1 



N f - Ni In 10 a A 



[exp lo (a A 0i + 6a) - exp lo (a A 0mm + 6 A )] ; (44) 



where, in any case, the maximum value relates to m ; n = 0. 

The fractional number of long-lived stars with normalized oxygen abun- 
dance, (pf < (p < max , can be inferred from Eq. (TIB"]) as: 

(MPs _ 1Qb , p" 

N f - Ni (Zoh 



x 



[l-^J^/^-a-PE^)^/^')] ; (45) 



or, in the linear limit (c — > 0), from Eq. fT23|) as: 

(SN) E 1 1 



N f - Ni In 10 a E 



[exp lo (a E 0max + 6 E ) - exp lo (a E 0/ + 6 E )] ; (46) 



where, in any case, the maximum value relates to max = ^0/(^0)0 via 
Eq. (18]) , which implies Z = 1 and may be replaced, with unrelevant numerical 
change in the results, by max —> +00. 
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Table 9: Upper values of the fractional star mass change due to the low- 
metallicity tail (A), the high-metallicity tail (E), both tails (A+E), under the 
assumption that the linear fit to the EDOD maintains unchanged outside the 
oxygen abundance range exhibited by the parent sample. Calculations were 
performed in the linear limit using values listed in Tables El and H] together 
with min = and max = Ao/(Z ) Q i- e - Z = 1 which, for practical purposes, 
is equivalent to max —> +00. See text for further details. 



(6sf) A /s f 




{Ss f ) E /s 


/ 


(5s f ) A+E /s f 


sample 


4.1925 E- 


-3 


2.8929 E- 


-3 


7.0854 E-3 


PMlln 


6.6414 E- 


-2 


1.1666 E- 


-2 


7.8079 E-2 


Ra07n 


1.8307 E- 


-1 


1.8189 E- 


-2 


2.0126 E-l 


PMllk 


3.1718 E- 


-2 


7.4427 E- 


-2 


1.0614 E-l 


Ra07k 



The fractional star mass change at the end of evolution, due to the pres- 
ence of (undetected) low-metallicity and high-metallicity tail in the EDOD, 
can be inferred under the assumption that star number is proportional to 
star mass, 5N/N = Ssf/sf. The result is: 

6^ (5N) A + (5N) E 

s f N f - Nl ' [ () 

where the values of the parameters can be deduced from the linear fits to the 
EDOD, as shown in Tables [3] and HI The results are listed in Table El with 
regard to upper values, in the linear limit. 

It can be seen the larger contribution arises from the low-metallicity tail 
for PMllk sample and from the high-metallicity tail for Ra07k sample. The 
reason is that, in both cases, the paucity of data makes more extended metal- 
licity tails, as shown in Fig.[2j The total change in fractional star mass does 
not exceed about 10% at most, with the exception of the (poorly populated) 
PMllk sample, where it raises to about 20%. Accordingly, the results of the 
current paper hold to an acceptable extent even if the linear trend shown by 
the EDOD is extended to undetected metallicity tails. 

In general, the disk is usually conceived as made of two main subsystems: 
the thick disk and the thin disk. Accordingly, the EDOD related to the disk 
depends on the thick to thin disk mass ratio, Mk/Mn, which is poorly known 
at present. Values already quoted in literature span a wide range, from some 
percent (e.g., Holmberg et al. 2007) to about unity (e.g., Fuhrmann 2008), or 
even indeterminate in the sense that no distinction can be made (e.g., Norris 
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1987; Ivezic et al. 2008; Bovy et al. 2012). Very low values of M K /M N could 
be biased unless corrected taking into consideration the different heights of 
the two subsystems on the Galactic plane at the Sun. 

The idea of a thick disk - thin disk collapse is not in contradiction with 
related specific angular momentum distribution, as shown in earlier attempts 
(Wyse and Gilmore 1992; Ibata and Gilmore 1995). If it is the case, in the 
light of the model, the thick and thin disk gas and star fractional mass are 
normalized to the same initial mass and the mass ratio of the thick to the 
thin disk at the end of evolution reads: 

m n (n f ) K + 0/) N + 0/) N 

where s/ and fif relate to the end of evolution and the indexes, K and N, 
denote the thick and the thin disk, respectively. From the results listed in 
TableEJ it can be inferred M K /M N = 0.08, 1.89, for PM11 and Ra07 samples, 
respectively. The former value would argue for a thick disk - thin disk col- 
lapse, but the PMllk sample is poorly populated to draw firm conclusions. 
On the other hand, the latter value would be against a thick disk - thin disk 
collapse, unless a massive thick disk (with respect to the thin disk) is recog- 
nized. Finally, if a single (thick + thin) disk population exists, the EDOD 
maintains close to its counterpart inferred both from the PMlln sample, due 
to the smallness of the PMllk sample, and from the Ra07n sample, as shown 
in an earlier attempt (Caimmi and Milanese 2009). 

The family of TDOD curves, considered for each stage of evolution, 
is characterized by two common points, (0j, ■?/>«), (0/,^/), which make the 
boundary of the related linear fit. As shown in Table [7] for the PMlln sam- 
ple, A stage, curves with increasing cut, Co, exhibit a monotonic trend for 
the remaining parameters, with the exception of the star mass fraction at the 
end of evolution, Sf, which attains a maximum. Accordingly, for a selected 
stage of evolution and fixed starting and ending point of the related linear 
fit to the EDOD, the star formation efficiency cannot exceed a threshold. In 
the case under consideration, the model predicts an amount not exceeding 
3.63 the initial mass of the thin disk was turned into stars during the earlier 
(A) stage of evolution. 

7 Conclusion 

The main results of the current paper may be summarized as follows. 

(1) The empirical differential oxygen abundance distribution (EDOD) in- 
ferred from two different samples for both the thin and the thick disk is 
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consistent with a linear trend within four, three or two regions (Fig. [2]). 

A family of theoretical differential oxygen abundance distribution 
(TDOD) curves, passing through the starting point, (0$,^), and the 
ending point, (0/, ■?/>/); of the linear fit to the EDOD, is defined within 
the framework of simple multistage closed (box+reservoir) (MCBR) 
models. In addition, regions where the EDOD exhibits a linear trend 
are related to stages of evolution characterized by different inflow or 
outflow rate and/or different oxygen abundance within the flowing gas. 
The extreme curves of the family on the (000) plane are related to 
Co — 0, (Co)oo) the last corresponding to k — > — oo, for inflowing gas 
and Co = (Co)ooj Ao/(Zo)f, the first corresponding to k — > +oo, for 
outflowing gas. 

For a family of TDOD curves related to an assigned stage of evolution, 
all the parameters show a monotonic trend passing from the lower to 
the upper curve, with the exception of the fractional stellar mass at 
the end of evolution, which attains a maximum and, in consequence, 
cannot exceed a threshold. In particular, the model predicts an amount 
not exceeding 3.63 the initial mass (Table [7]) of the thin disk is turned 
into stars during the earlier (A) stage of evolution, in connection with 
the EDOD inferred from the PMlln sample. 

The special case of steady inflow regime, where the TDOD reduces 
to a horizontal line, is consistent with the results of hydrodynamical 
simulations, where quasi equilibrium is attained between inflowing gas, 
outflowing gas and gas lost via star formation (e.g., Finlator and Dave 
2008; Dave et al. 2011a,b; 2012). 

If the linear trend inferred for the EDOD extends towards undetected 
low and high oxygen abundance, then the fractional star mass change 
does not exceed about 10% for well populated samples (PMlln, Ralln, 
Rallk) and about 20% for poorly populated samples (PMllk). 

Under the assumption of a thick disk - thin disk collapse, model pre- 
dictions yield a mass ratio, M K /M N = 0.08, 1.89, with regard to PM11 
and Ra07 samples, respectively. The latter alternative is in contradic- 
tion with current observations, which seem to exclude the presence of a 
massive thick disk, while the former remains available on this respect. 
If thick disk and thin disk stars belong to a single population (e.g., Bovy 
et al. 2012), no prediction can be made for or against a thick disk - thin 
disk collapse and the TDOD is expected to be slightly different from 
its counterpart related to the thin disk (Caimmi and Milanese 2009). 
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The main uncertainties on the EDOD used in the current paper are related to 
(i) poorly populated samples (PMllk); (ii) biased samples (all) towards low 
metallicities, [Fe/H] < —1.0 or [O/H] < —0.5, which are expected to affect 
mainly the thick disk, where the overall metallicity or oxygen abundance are 
lower than in the thin disk, even if the low-metallicity tail is neglected. On the 
other hand, the prediction of early strong inflow implies the main features 
of thick and thin disk evolution are captured by the EDOD inferred from 
PMlln, Ra07n, Ra07k samples, even if the low-met allitity tail is neglected. 
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Figure 1: The empirical, differential oxygen abundance distribution (EDOD) 
related to PM11 (upper panels) and Ra07 (lower panels) samples. The un- 
certainty of the distribution is determined from Poisson errors, which tend 
to negative infinite for data related to a single star, including the one outside 
the bottom left panel. The vertical dotted lines mark the boundary between 
adjacent regions, as determined from the intersection of related linear fits. 
For further details refer to the text. 
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Figure 2: Regression lines to the empirical differential oxygen abundance 
distribution (EDOD) plotted in Fig.[H with regard to the regions (from the 
left to the right): A, F, C, E. If adjacent regions exhibit a similar linear 
trend, a single fit is performed. Data related to a single star are not used 
for the fitting procedure unless only two points remain. Other captions as in 
Fig.[TJ For further details refer to the text. 
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Figure 3: The parameter space on the (Ok(o) plane. The locus, c = 0, 
is represended by the equilateral hyperbola lying on the second and fourth 
quadrant, where c > or c < for points placed between the branches 
of the hyperbola or outside either branch, respectively. The locus, /i/ = 
0, is represended by the equilateral hyperbola lying on the first and third 
quadrant, where /i/ > or fif < for points placed between the branches of 
the hyperbola or outside either branch, respectively. The vertical axis, k = 0, 
and the horizontal axis, Co — 1? are asymptotes for both hyperbolas. The 
vertical axis, k = —1, represents the steady inflow regime. The parameter 
space is restricted to the region bounded by the branches of the hyperbola, 
/if — 0, and the horizontal axis, Co = 0, which ensures /if > and Co > 0. 
For further details refer to the text. 
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Figure 4: Comparison between empirical (EDOD), plotted in Fig.HJ and the- 
oretical (TDOD) differential oxygen abundance distribution defined by the 
family of curves including the regression line and passing through the starting 
point, ((pi, ipi), and the ending point, (</>/, ipf), where ip u = a<p u + &n, u = i, f, 
&n = b + logC*N, for the stage under consideration. Curves correspond (from 
bottom to top) to cut values, Co, between and (Co)oo or (Co)oo and 1.2 for 
negative and positive flow values, k, respectively. With regard to the top 
left panel, the lower curve fitting the data relates to Co — 0.35. In general, 
for varying cut, (Co)min < Co < (Co)max, the TDOD lies within the region 
bounded by the curves, Co = (Co)min and Co = (Co)max- The full straight 
line on the top of each panel has slope, ao = — (1/ In 10)[(Zo)o/p] = —0.42, 
and relates to null flow, k = 0, thus separating the inflow regime, a > a , 
k < 0, from the outflow regime, a < a Q , k, > 0. A horizontal straight line 
corresponds to the steady inflow regime, a = 0, k = — 1. 



31 



